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Abstract 

The two-band Hubbard model is used to analyze a possibility of a non-uniform charge distribution 
in a strongly correlated electron system with two types of charge carriers. It is demonstrated 
that in the limit of strong on-site Coulomb repulsion, such a system has a tendency to phase 
separation into the regions with different charge densities even in the absence of magnetic or any 
other ordering. This tendency is especially pronounced if the ratio of the bandwidths is large 
enough. The characteristic size of inhomogeneities is estimated accounting for the surface energy 
and the electrostatic energy related to the charge imbalance. 

PACS numbers: 75.47.Lx, 64.75. +g, 75.30.-ni, 71.30.+h 
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I. INTRODUCTION 



The phase separation is commonly considered as an inherent property of strongly corre- 
lated electron systems^. Usually this phenomenon is treated as a result of a coexistence and 
competition of different kinds of ordering (magnetic, charge, orbital)^'-^^. The most widely 
discussed type of the phase separation is a formation of nanoscale inhomogeneities such as 
ferromagnetic metallic droplets in an insulating antiferromagnetic material arising due to 
the self-trapping of charge carriers^. Such type of the phase separation is characteristic of 
the doped manganites. The another type of the nanoscale inhomogeneity is the modulation 
of the electron density due to antiferromagnetic correlations, which is considered as possible 
mechanism of the phase separation observed in superconducting cuprates^i^. 

Nevertheless, the phase separation can manifest itself even without some specific type 
of ordering, e.g., if the system contains different types of charge carriers. The simplest 
illustration of such a behavior gives the Falicov-Kimball model^, which is often used as a 
toy model for heavy-fermion compounds. This model describes the system of itinerant and 
localized electrons with a strong on-site Coulomb repulsion. The numerical simulations of 
this model demonstrated an inhomogeneous charge density distribution at some relation 
between the itinerant electron bandwidth and the distance between the localized level and 
the bottom of the band^. The competition between the metallicity and localization in a 
similar system with magnetic interactions was studied in Refs. lolliol with an emphasis on 
the phase diagram of magnetic oxides such as manganites. The system with a band and 
localized level is a limiting case of much more common situation of two bands with different 
width. 

In this paper, we use a two-band Hubbard model for the description of a strongly cor- 
related electron system with two types of charge carriers. We demonstrate that the phase 
separation in this system arises even without any ordering if the ratio of the bandwidths 
is large enough. In Section II, we write out the Hamiltonian of the model. In Section III, 
we study the electron structure of a homogeneous state. In Section IV, we analyze the pos- 
sibility of the phase separation and estimate the size of inhomogeneity accounting for the 
long-range Coulomb interaction and surface energy. In Section V, we discuss the obtained 
results. 
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II. THE MODEL 

Let us consider a strongly correlated electron system with two bands a and h of different 
width. Let the first band, a, be wider than the second one, h. Such a system could be 
described by the following Hubbard Hamiltonian 

1 U' 

(ij}a,(T ic iQ,(T ia,(T ia^cjcr' 

Here, O;^^ and the creation and annihilation operators for electrons in the bands 

a = {a, b} at site i with spin projection a, and riiaa = alaa^^iaa- The symbol (...) denotes the 
summation over nearest-neighbor sites. The first term in Eq. ([T]) corresponds to the kinetic 
energy of the conduction electrons in bands a and b with the hopping integrals ta > ti,. In 
our model, we ignore the interband hopping. The second term describes the shift e of the 
center of band b with respect to the center of band a. The last two terms correspond to 
the on-site Coulomb repulsion of two electrons either in the same state (with the Coulomb 
energy f/") or in the different states {U'). The bar above a or a denotes not a or not a, 
respectively. The assumption of the strong electron correlations means that the Coulomb 
interaction is large, that is, f/", U' ^ t", e. The total number n of electrons per site is a 
sum of electrons in the a and b states, n = Ua + Ub, and fi is the chemical potential. Below 
for definiteness sake, we consider the case n < 1. 

III. HOMOGENEOUS STATE 

The homogeneous state of the model formulated above can be analyzed by standard 
methods at arbitrary band filling n. Let us introduce a one-particle Green function 

Gaa{} - jo, t - to) = -i{fajaa{t)a]^^^{to)), (2) 

where T is the time ordering operator. The equation of motion for the one-particle Green 
function with Hamiltonian ([T]) can be written as 

(^^^ + /i + G^^a - jo, t - to) = 5jjo5(t - to) - t" J^G.^a - jo + A, t - to) 

+ U'^Qaa,aa (j " jo, t - to) + U' ^ Gaa.aa' (j " jo, t - to) , (3) 

a' 
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where e"" = for a = a and = e for a = b, the summation in the second term in the right- 
hand side of Eq. ([3]) is performed over sites nearest to j, and A are the vectors connecting 
the site j with its nearest neighbors. Equation ([2]) includes the two-particle Green functions 
of the form 

Qaa,f3a'{^ " jo, ^ " ^o) = Ojc^^ (t)^j/3a' (t)^]^^.^ (^o) ) • (4) 

Then, we should write the equations of motion for these functions, which will include the 
next order Green functions, etc. To cut such infinite chain of equations, we shall use here 
the following procedure. 

In the limit of strong Coulomb repulsion, the presence of two electrons at the same 
site is unfavorable, and the two-particle Green function Eq. (jl]) is of the order of 1/f/, 
where U ~ Ua,U'. The equation of motion for Qaa,i3cT' includes the three-particle terms 
coming from the commutator of ajaa{t) with the U terms of Hamiltonian ([1]) in the form 
{Tajaa{t)n^i3a' (t)nj^^'i {t)aj^^^(to)) . It is easy to see that these terms are of the order of 
and in our approximation, we neglect them. In the equations of motion for the two-particle 
Green functions, we make the decoupling corresponding to the Hubbard I approximation^i. 
That is, in term coming from the commutator of aj^ait) with the kinetic-energy terms 
of Hamiltonian ([1]), we make the following replacement {Taj+/^aa(t)nj[^„i{t)aj^^^^(to)) — >■ 
(^j/3(T')(^'^j+AQ:cr(i)ctjf,Q,o-(io))- The analogous decoupling in terms coming from the com- 
mutator of rijaait) with the same kinetic-energy operator yields zero^^^^^. As a result, the 
equations for the two-particle Green functions can be written as 



d 

H + e°' -U°' ] ^a<T,as(j - jo, ^ - to) 



hjii - to) - - jo + A, t - to) 



(5) 



+ A* + ~ ^' ) Gaa,aM - jo, ^ " ^o) 



Tin 



^iio^it - to) -^Yl ^--(j - jo + A, t - to) 

A 



(6) 



where = {^jaa) is the average number of electron per site in the state (a, a). 

Eqs. ([3]), and (jS]) are the closed system for one- and two-particle Green functions. 
This system can be solved in a conventional manner— passing from the time-space to the 
frequency-momentum representation. Eliminating two-particle Green functions we can find 
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the explicit expression for single-particle Green functions, which have poles corresponding 
to the Hubbard sub-bands for each band, a or h. The splitting between these sub-bands is 
determined by on-site Coulomb repulsion of electrons in the same states. The on-site 
Coulomb repulsion U' of electrons in different states gives rise to the correlation between the 



a and h bands (analogous to the case of localized and itinerant electrons discussed 



fillings o 
in Ref.lgllOl). 

If the total number of the electrons per site does not exceed unity, n < 1, the upper 
Hubbard sub-bands are empty, and we can proceed to the limit f/' ^ oo. In this 
case, the one-particle Green function Gaa in the frequency-momentum representation can 
be written as 

^-(k' ^) = ^ ^ /"'^ TriA ' 

where Wa = zt"', z is the number of nearest neighbors, 

fl'a<7 = 1 - ^ naa' " , (8) 



and 



is the spectral function depending on the lattice symmetry. In the case of simple cubic 
lattice, we have 

((k) = -- [cos{k^d) + cos{k^d) + cos{k^d)] , (9) 
3 

where d is the lattice constant. 

In the main approximation in 1/f/ considered here, the magnetic ordering does not ap- 
pear. To study the possible types of magnetic ordering in our model, it is necessary to take 
into account the higher order terms in 1/U. So, below we assume that 

"-aT = ^ai = rial 2. (10) 

From the expression for Green function ([7]), it follows that the filling of each lower sub-band 
is equal to Qa] + Qai = ^g^. Using ([HD and (ITUi) . we have 

5f„ = 1 - - y . (11) 

Using the relation for the density of states Pa{E) = — vr^^Im J Ga(k,E + i0)(i^k/(27r)^, 
we get the following expression for the number of electrons 

Ha = ^gcJiQ (12) 

V 9aWcx J 




Figure 1: (Color online) Ua (solid line) and (dashed line) vs the total number of charge carriers 
n; Wb/wa = 0.2, e/wa = 0.12. Vertical arrows show the concentrations of and Ub in the 
inhomogeneous state (see the text below). 

where 

no(/i') = J dE'po{E'), (13) 
-1 

and 

Po{E') = / ^(^' - C(k)) (14) 

is the density of states for free electrons (with the energy normalized by unity, \E\ < 1). 
The chemical potential /i in Eq. (fT2l) can be found from the equality n = na + Ub. 

Let us consider the case when the energy difference e between centers of a and b bands 
is not too large, that is, of the order of the width of b band, Wb- In this case, there exist 
only a electrons at low doping until the chemical potential reaches the bottom of the b band 
— e — Wb at some concentration Uc- At n > ric, the b electrons appear in the system, and the 
effective width of a band, Wa = 2waga{na,nb), starts to decrease. The plots of n^, Ub, and 
the effective bandwidth are shown in Fig. [1] and Fig. [2] as functions of n. In all calculations, 
we use the spectrum ^(k) in the form ([9]). 



6 



-1.0 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



n 



Figure 2: (Color online) Effective bandwidths Wa = 2waga vs the total number of charge carriers 
n. The dashed curve is the chemical potential fi. The values of the parameters are Wb/wa = 0.2, 
e/wa = 0.12. Vertical arrows show the concentrations of Ua and rih in the inhomogeneous state. 

The energy of the system in homogeneous state, -Ehom, is the sum of electron energies in 
all filled bands. Using the the density of states for free electrons Eq. (|T^ we can write -Ehom 
in the following form 




(15) 



where 




(16) 



-1 



The dependence of -Ehom(^) is shown in Fig. [3] by the solid line. 



IV. PHASE SEPARATION 



A. 



General consideration 



In this section, we analyze the possibility of the phase separation in the system. As 
one can see from Fig. [3l the energy for the homogeneous state, E^omin), has two minima 



at different values of the charge carrier density. In this situation, it is favorable for a 
system to form two phases with different electron concentrations. Moreover, the existence 
of two minima is not a necessary condition for the formation of an inhomogeneous state and 
this phenomenon could be observed under more general conditions^i^^. However, the phase 
separation could be hindered by increase of the energy due to surface effects and a charge 
redistribution. 

At first, we do not take into account the charge redistribution and surface terms in the 
total energy. In this way, we determine the optimum content of each phase and the corre- 
sponding charge densities but we can not find the characteristic size of the inhomogeneities. 

We consider two phases, I (low carrier density) and II (high carrier density), with the 
number of electrons per site ni and n2, respectively. A fraction p of the system volume is 
occupied by the phase I and 1 — p is a fraction of the phase II. We seek a minimum of the 
system energy 

E°^{ni,n2) = pEi^omini) + (1 - p)-Ehom(^2) (17) 
under the condition of the charge carrier conservation 

n = prii + (1 — p)n2 ■ (18) 

The results of calculations of the system energy in the phase separated state are shown in 
Fig. [3]for Wb/wa = 0.2 and e/wa = 0.12 by the dashed curve. From this figure we see, that 
the phase separation exists in the range of n where both types of charge carriers coexist 
in the homogeneous state. The numerical analysis shows that the concentrations of the 
charge carriers in each phase, ni and n2, vary slowly with n, remaining close to certain 
optimal values for each phase: ni ~ 0.5 for the state with low carrier density, whereas 
n2 ~ n;, ~ 1 for the state with high carrier density. The phase II can be considered as a Mott- 
Hubbard insulator since the corresponding lower Hubbard sub-band is almost completely 
filled. If n increases from zero to one, the phase separation may be favorable when the n 
achieves the value corresponding to the energy minimum for the homogeneous state. At this 
value of n the content of the phase I p{n) starts to decrease from p{n) = 1. With the further 
increase of n, p{n) tends to zero at ^ Therefore, we can conclude that the system may 
separate into metallic and insulating phases in a certain parameter range. An indication to 
the phase separation is a negative curvature of the Ehom{n) curve at the right side from the 
energy minimum, see Fig. [3l 
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The above discussion demonstrates that the width and filhng of one band depends on the 
width and filhng of other bands. The phase separation gives the possibihty to attain the 
minimum free energy by an optimum filhng of the electron bands. The phase separation can 
be favorable only if the bands are appreciably different. In Fig. |3]we plot the energy of the 
homogeneous state versus n at different values of ratio Wh/wa- We see that d'^E /dn^ < 
in a wide range of n ioi Wb/wa ^ ^, which indicates the possibility of the phase separation 
(see also Fig. [3]). The phase separation becomes unfavorable only if Wb/wa < 0.4. 



B. Characteristic size of inhomogeneities 



The phase separation leads to redistribution of charge carriers {rii ^ n2). Therefore, 
we should take into account the electrostatic contribution to the total energy of the phase- 
separated state. This contribution depends on the shape of inhomogeneities. Following 



Refs. 
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we consider a spherical geometry of the inhomogeneous state, namely, at 
p < 0.5, the sample is modelled as an aggregate of spheres of phase I embedded into the 
matrix of phase II and vice versa for p > 0.5. The electrostatic, Eq, energy is calculated in 
the Wigner-Seitz approximation, that is, we consider a set of spherical unit cells with zero 
total charge, where the spherical core of one phase is surrounded by a shell of another phase. 



Following the approach described in Refs. 
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Er 



5ed 



{ni - n2)' 



we find at p < 0.5 



(19) 



where e is the average permittivity and Rg is the radius of the droplet of the phase I. In the 
case p > 0.5, we should replace ni ^ n2 and p ^ 1 — p. 

Another contribution to the total energy depending of the size of inhomogeneities is 
related to the surface between two phases. It comes from the size quantization and depends 
on the electron densities in both phases. The case when one of the densities is zero was 



considered in Ref. 



id . The generalization of this approach for non-zero densities is presented 



in Appendix, where surface energy a{ni,n2) is calculated using the perturbative approach 



proposed in Ref. 



13l . The corresponding contribution, Es-, to the total energy is proportional 



to the surface area between phases I and II. hi p < 0.5, it can be written in the form 

3ci , 



E. 



p—a[ni,n2) 

/to 



(20) 
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Figure 3: (Color online) The energy of the system vs doping level n. Solid curve corresponds to the 
homogeneous state, whereas dashed curve is the energy of phase-separated state without taking 
into account electrostatic and surface contributions to the total energy. Dot curves are the energies 
of inhomogeneous state, Eqs. ([22]) and ([23|) . at Vo/wa = 0.1, 0.05, 0.01 from top to bottom (see 
text below). Here, Wh/wa = 0.2 and e/wa = 0.12. Vertical arrows show the concentrations of 
and rib in the inhomogeneous state. 

In the case p > 0.5, we should replace p ^ 1 — p. 

Minimization of the sum i?cs = Eq + with respect to Rs allows us to calculate this 
value. In doing so, we get at p < 0.5 

15a{ni,n2) (21) 



(22) 



'4nVo{n2 - ni)^ (2 - 3p^/^ + p) 
The total energy of the inhomogeneous state then reads 



Eps{ni, ris) = pEhom{ni) + (1 - p)£'hom(^2) + £^cs(^i, ^2) 



where 



97r 



Ecs = s( Vo — {n2 - ni)V(ni,n2) ) p {2 - 3p'/'' + p) [ 



1/3 



1/3 



(23) 



Vb = (? jed. The energy i?ps calculated by minimization of Eq. (!22l) with respect to n\ and 
77-2 at different values of Vq are shown in Fig. [3] by dot-dashed curves. We see, that the 
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n 

Figure 4: (Color online) The energy of homogeneous state vs doping level n at Wh/wa = 0.1 (solid 
line), Wb/wa = 0.25 (dashed line), and Wb/wa = 0.4 (dot-dashed line). The parameter Sb/wa = 0.1 
for all cases. The phase separation is favorable for solid and dashed curves in the range of doping 
0.45 ^ n < 1, where d'^Ejdr? < 0, whereas for dot-dashed curves, only the homogeneous state 
exists. In the inset, the maximum energy gain due to the formation of the inhomogeneous state 
(Vo = 0, Sb/wa = 0.1) as a function of the ratio Wb/wa is shown. For Wb/wa ^ 0.38 the phase 
separation becomes unfavorable in energy. 

Coulomb repulsion in the system with an inhomogeneous charge distribution reduces the 
range of n, in which the phase separation is favorable. In Fig. O we plot the characteristic 
radius of inhomogeneities, i?^, function of n. 

V. DISCUSSION 

Thus, the phase separation can be favorable for the system of the strongly correlated 
electrons even in the absence of any specific ordering. We demonstrated that the state with 
inhomogeneous charge distribution can arise if there exist two types of the charge carriers 
with different bandwidth. The electron correlations due to on-site Coulomb repulsion lead 
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Figure 5: (Color online) The radius of droplets Rs vs doping level n at Wb/wa = 0.2, e/wa = 0.12, 
and Vo/wa = 0.02. 

to the dependence of the bandwidth for one type of electrons on the band filling for another 
type of electrons. As a result, the dependence of energy on the total number of the charge 
carriers becomes non-monotonic. The competition between kinetic and correlation energies 
triggers the formation of an inhomogeneous ground state. It is particularly evident if the 
energy of the system as a function of electron density has two minima. In this case, it could 
be favorable for the system to separate into two states with electron densities close to these 
minima rather than to form a homogeneous state with an intermediate density. Such a 
situation is illustrated in Fig. [31 

It is clear that the phase separation can occur only if the bandwidths corresponding to 
two types of the charge carriers are sufficiently different, that is, the ratio Wb/wa of the 
widths of narrow and wide band should be rather small. The second condition is that the 
narrow and wide bands should not be widely separated from each other, that is, the ratio 
e/wa of the distance between the band centers and the width of the wider band should be 
less than unity. Naturally, the long-range electrostatic forces prevent nonuniform charge 
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distribution and the condition Vo/wa ^ I should be met. As it can be seen from Fig. HI the 
phase separation can be favorable even if Wb/wa < 0.3 — 0.4. 
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Appendix: SURFACE ENERGY 

The electrons in the phase-separated state are confined within a restricted volume Vg. 
This gives rise to the change in the density of states in both phases. At small ratio A = 
Ssd/Vs, where Sg is the surface area of the sample, the density of states for free electrons 
can be written as^>i^ 



A 



A 



1+y iPo(i5')- j(pr^(^'+i/3)+pr^(i5^'-i/3)) . (A.i) 



where Pq is given by Eq. (fT4l) . Here p'^^^ is the density of states in two dimensions. Using 
this expression instead of f|T4l) . and expanding the Eqs. f|T2l) . f|T3l) . f|T5l) . and f|T6l) in a series 
in powers of A up to the first order, we derive formula (l20l) with the correction for the size 
quantization (j{n) in the form 



+ 



2p 



(1) 



Wr 



H'a "'a 



en 



(1) 

b ' 



(A.2) 



where 



„(0) I .a 
'(0) _ P + ^ 



P 



(0) 



n 



(1) ^ 

~ 2 



'2no{pT)-nr\pT' + l)-nr\pT'-l) 



1 + Pa^ Po{pa^) - noipa''') 

1 

' 12 



,'(0)^ 



feo(p') = ^£o(p') - 2 



(A.3) 



(A.4) 



n?'\p' + l)-nr{p'-l) 



(A.5) 
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correction to the chemical potential /i*^^-* is found from the condition Y2a ~ 0' ^"^^ 
superscript (0) denotes the unperturbed value of corresponding quantity. The functions 
n^Q^\jj') and e''^^\jj') in these expressions are determined by Eqs. f|T3l) and f|T6l) . respec- 
tively, where one should change po to Po^^^ The surface energy a{ni,n2) per unit area 
between phases I and II is the sum a{ni) + a{n2). 
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